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The measured hourly daily data of total, direct and diffuse solar irradiation incident on a horizontal and an 
inclined surface for Cairo, Egypt (Lat. 30°05’N and Long. 31°15’E), during the period (1990-2010) are analyzed. 
The regression equations between (G/G,) and meteorological variables along with the values of MBE, RMSE, 
MPE, R? and the t-test statistics are summarized in this research. The values of correlation coefficients (R?) are 
higher than 0.95 and the values of the RMSE are found in the range 3.13-6.34, thus indicating a good 
agreement between measured and calculated values of the total solar radiation (G). The models of Eqs. (10), 
(11) and (14) have well estimated the total solar irradiation in the selected location during the time period in 
the present study. For all models, the absolute values of the MPE indicate very good agreement between 
measured and calculated values of the diffuse solar fraction (Gg/G) or the diffuse solar transmittance (Gg/G,) 
and clearness index K, relative number of sunshine hours (S/S,) and their combination. The models of Hay 
(Ha), Skartveit and Olseth (SO) and Perez et al. (P9) give the most accurate predictions for the south-facing 
surface, and Hay (Ha) and Perez et al. (P9) models performs better as estimated for the west-facing surface. 
© 2013 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Solar irradiation is measured at many locations around the 
globe. Unfortunately, these locations are mainly concentrated in 
developed countries, and are scarce within the developing world; 
solar irradiation data are fundamental inputs for solar energy 
applications such as photovoltaics, solar thermal systems, and 
passive solar design. Daily data on horizontal and tilted surface are 
required for day-to-day performance monitoring in applications 
such as solar gains from vertical glazing, day lighting and agricul- 
tural processes; these data are also useful to engineers in design- 
ing of various solar energy conversion devices and accurate 
modeling of the impact of solar gains through glazing which is 
imperative especially when simulating the thermal behavior of 
these buildings. Empirical validations of solar gain models are 
therefore an important and necessary endeavor to provide con- 
fidence to developers and modelers that their respective algo- 
rithms simulate reality [1-3]. 

Solar energy consists of two parts: extraterrestrial solar energy 
which is above the atmosphere and global solar energy which is 
under the atmosphere. The measured solar energy values can be 
used for developing solar energy models which describe the 
mathematical relations between the solar energy and the meteor- 
ological variables such as ambient temperature, humidity and 
sunshine ratio. These models can be later used to predict the 
direct and diffuse solar energy using historical metrological data at 
sites where there is no solar energy measuring device installed. 
Solar energy of daylight utilization for any site is dependent upon 
the quality of the available flux. Obviously, the flux impinging 
upon any arbitrary surface undergoes monthly as well as diurnal 
variations. The measurements of solar energy received from the 
sun, on horizontal as well as on sloped surface, are an expensive 
affair. As such, few locations in the world have reliable, long-term 
measured irradiation data sets. Daylight records are even scarcer. 
Most radiation data are given as the energy received on a 
horizontal surface. Since only very few applications use this 
configuration, there is a genuine need for insolation estimation 
to be carried out for sloped surface of any given aspect; the 
accuracy of these models varies from 40% to 50% for abbreviated 
techniques to limits set out by the accuracy of the measuring 
equipment for modern sophisticated models [4,5]. The total solar 
radiation on a horizontal surface is called global irradiance and is 
the sum of incident diffuse radiation plus the direct normal 
irradiance projected onto the horizontal surface. If the surface is 
tilted with respect to the horizontal, the total irradiance is the 
incident diffuse radiation plus the direct normal irradiance pro- 
jected onto the tilted surface plus ground reflected irradiance that 
is incident on the tilted surface [6-9]. 

Many empirical models have been used to estimate solar 
radiation, utilizing the available data on meteorological [10,11], 
geographical and climatological parameters. Among these para- 
meters, sunshine duration [12-15], air temperature [16], latitude 
and longitude [17], precipitation [18], relative humidity [16,19,20], 
wind speed and cloudiness [21-26] were used. Of these, the most 
commonly used parameter for estimating total solar radiation is 
sunshine hours. In this respect, the modified version of Angstrom 
equation, among various correlations, has been widely used to 
estimate the total solar irradiation on horizontal surface [27-35]. 

In most of the solar energy applications, inclined surfaces at 
different angles are widely employed. The solar irradiance on a 
horizontal surface has been measured in many meteorological 
stations around the world, but there are only a few stations that 
measure the solar component on inclined surfaces [36-38]. There 
are a number of models available to estimate solar irradiation on 
inclined surface from corresponding horizontal data. This requires, 
in general, the availability of detailed information on the 


magnitude of diffuse and direct horizontal irradiance. A number 
of diffuse fraction models are available as documented in [39-43]. 
These models are usually expressed in terms of polynomial 
functions relating the diffuse fraction to the clearness index. An 
inclined surface solar irradiation model developed by Olmo et al. 
[44] requires only the horizontal surface solar irradiation, with 
incidence and solar zenith angles as input parameters. 

Many solar energy models have been presented in the litera- 
ture using mathematical linear [45-53] and nonlinear functions 
[54-60], artificial neural network [61-71] and fuzzy logic [72-74]. 
An important aspect in modeling solar energy is the accuracy of 
the developed model which is evaluated using statistical errors 
such as the mean absolute percentage error (MAPE), mean bias 
error (MBE) and root mean square error (RMSE). The MAPE is an 
indicator of accuracy in which it expresses the difference between 
real and predicted values to the real value. The calculated MAPE is 
summed for every fitted or forecasted point in time and divided 
again by the number of fitted points; n. MBE is an indicator for the 
average deviation of the predicted values from the measured data. 
A positive MBE value indicates the amount of overestimation in 
the predicted total solar energy and vice versa. On the other hand, 
RMSE provides information on the short-term performance of the 
model and is a measure of the variation of the predicted values 
around the measured data. RMSE also shows the efficiency of the 
developed model in predicting future individual values [3,71]. A 
large positive RMSE implies a big deviation in the predicted value 
from the measured value. Part of this study was used before to 
evaluate the statistical comparison models of solar energy on 
horizontal and inclined surfaces [3]; in the present work, the 
average hourly daily data of solar irradiation models on horizontal 
surface and various inclined surface in the selected site will be 
studied, using empirical models which were selected to estimate 
the solar radiation on horizontal and inclined surfaces. 


2. Instrumentations and climate site 


In the present work, the global, direct and diffuse solar 
radiation incident on a horizontal surface at Cairo, Egypt (Lat. 
30°05'N and Long. 31°15’E), during the time period from January 
1990 to December 2010 is used. The radiation data of the 
corresponding periods are obtained from the Egyptian Meteor- 
ological Authority. The data sets used consist of mean hourly and 
daily values of global and diffuse solar irradiations on a horizontal 
plane. Total solar radiation was measured using Eppley high- 
precision pyranometer responsive at 300-3000 nm, while another 
precision pyranometer equipped with a special shading device, 
SBS model, was used to measure diffuse irradiation. The shadow 
band stand is constructed of anodized aluminum, which weighs 
approximately 24 lb and uses a 300 band of approximately 2500 
diameter to shade the pyranometer. Because the shadow band 
screens the sensor from a portion of the incident diffuse radiation 
coming in from the sky, a correction was made to the measure- 
ments following Batlles et al. [3,75,76]. Total solar radiation data 
were recorded by the Eppley Precision Spectral Pyranometer (PSP) 
at all stations. The accuracy of these pyranometers corresponds to 
the first class according to the World Meteorological Organization 
classification [77]. These instruments are calibrated each year 
against a reference instrument traceable to the World Radiometric 
Reference (WRR) maintained at Davos, Switzerland [78,79]. 
According to the calibration certificate of the manufacturers, 
sensitivity is approximately 9 „V/W/m?, temperature dependence 
is +1% over ambient temperature range —20 to +49 °C, linearity 
is + 0.5% from 2800 W/m”, and cosine is + 1% from normalization 
0-70° zenith angle and + 3% for 70-80° zenith angle. The absolute 
accuracy of calibration is + 3-4%. 
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There are few peculiarities in the climate of Egypt that may 
affect the characteristics of the sough relation, which we shall 
outline here. Cloud characteristics and air temperature change 
from one season to another as follows: winter is the season of 
different cloud types, which are normally opaque to the direct 
beam, and has minimum temperature, in addition to low turbidity 
of the atmosphere. Spring is characterized by the passage of small 
and shallow thermal depressions, inducing what is called Khama- 
sin weather. Vertical visibility is deteriorated progressively with 
increasing dust content in the lower layers. After the passage of 
these depressions, clouds form and horizontal visibility is reduced 
much in the atmosphere. In summer, high temperature, high 
transparency and semi-transparent clouds prevail, but despite 
their existence, the sky is ‘dirty’ most of the time, due to deep 
layers of fine dust particles associated with continental tropical air. 
The dust content falls markedly when Mediterranean air arrives, 
associated with fine weather cumulus. In autumn, the atmosphere 
is fairly transparent on average. Morning mists and low clouds 
dissipate after sunrise [3,80]. 


3. Solar radiation basic 


3.1. Computing the semi-hourly total extraterrestrial solar radiation 
Gon 


Solar radiation incident outside the earth's atmosphere is called 
extraterrestrial solar radiation. On average the extraterrestrial 
irradiance is 1367 W/m? (solar constant). The extraterrestrial 
radiation Gon is given as follows [3,81]: 


Gon = (24/7) x Isc 
xE [cos pcos sin w+ (t@/180)sin gsin ô] (1) 


where E, is the correction factor of the Earth's orbit and œ is the 
sunrise/sunset hour angle given respectively by 


E» = 1 + 0.033cos(2zdn/365), (2) 


w= cos 7! (— tan o tan ô) (3) 


where ¢ is the latitude and the solar declination angle of the sun 
(6) is the angle between a plane perpendicular to a line between 
the earth and the sun and the earth's axis, which is given in 
degrees according to Spencer [82] as 


5 = (0.006918—0.399912cosr + 0.070257sinr—0.006758cos2r 
+ 0.000907sin27—0.002697cos3r + 0.00148sin3/)(180/z) (4 


where I is the day angle in radiance, and it is represented by 
T = 2n(dy—1)/365 (5) 


where dn is the day of the year. 


3.2. Zenith azimuthal and hour angles 


To describe the sun's path across the sky one needs to know the 
angle of the sun relative to a line perpendicular to the earth's 
surface which is called the zenith angle (@) and the angle of the 
sun's position relative to the north-south axis is the azimuthal 
angle (a). The hour angle (œ) is easier to use than the azimuthal 
angle because the hour angle is measured in the plane of the 
“apparent” orbit of the sun as it moves across the sky (Fig. 1) [3,6]. 
With the above information, one can calculate the cosine of the 
zenith angle as follows: 


cos (Z) = cos (g) cos (6) cos (œw) + sin osin 6 (6) 


Zenith 


Horizontal 
Plane 


Fig. 1. Zenith, azimuthal, and hour angles. 


4. Modeling of solar energy 


The commonly used solar energy models developed in the past 
are based on linear and nonlinear models [83]. These models give 
a correlation between solar energy on a horizontal surface and 
some meteorological variables such as shine hours s, ambient 
temperature T, cloud cover Cw, relative humidity Rp, and maximum 
Tmax. and Tmin. ambient temperatures. The linear models use 
simple linear function while the nonlinear models are polynomial 
functions of the third or fourth degree. 


4.1. Global solar energy models on horizontal surface 


Linear and nonlinear models can be used to calculate the total 
solar energy in terms of sunshine hours. A commonly used linear 
model for this purpose which defines the total solar energy in 
terms of the extraterrestrial solar energy is given [3,5,84] as 
follows: 


G/Go =a + b(S/So) (7) 
G/Go =a + b(S/So) + C(S/So)* (8) 
G/Go =a + b(S/So) + cT (9) 
G/Go =a + b(S/So) + CRy (10) 
G/Go =a + bT + cRy (11) 
G/Go =a + D(Tmax.-Tmin.) + CCw (12) 
G/Go = a + D(Tmax.-Tmin.)”” + CCw (13) 
G/Go =a + b(S/So) + CCw 14) 
G/Go =a + b(S/So)® (15) 


where a, b and c are the empirical constants and sọ is the 
maximum possible monthly average daily sunshine duration or 
the day length. 


4.2. Computing beam and diffuse radiation on horizontal surface 


The empirical correlations for calculating the average hourly 
daily diffuse radiation incident on a horizontal surface, the diffuse 
fraction (Gq/G) and diffuse transmittance (Ga/Go) were correlated 
to first, second and third order correlations of the clearness index 
(K) and the relative number of sunshine hours (S/S,) [3,57]. It is 
found that the second and third order correlations do not improve 
the accuracy of estimation of (Ga). Therefore, the following 
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correlations have been obtained for Cairo: 


Ga/G = 5.817-6.517K;, R? = 0.972 (16) 
Ga/G = 8.342-6.455(S/So), R? = 0.958 (17) 
Ga/Go = 3.815-5.319(S/So), R? = 0.932 (18) 
Ga/Go = 4.912-6.894K;, R? = 0.985 (19) 


Furthermore, (Gqa/G) and (Ga/Go) were correlated to the first 
and second order correlations of the (K+) and (S/S,) combination. 
Also, it has been found that the second order correlations between 
(Gg/G) or (Ga/Go) and (K+) and (S/S,) combination do not improve 
the accuracy of estimation of (Ga). The following correlations were 
found to fit the measured data of (Gy): 


Ga/G = 6.314-5.131K; + 0.136(S/So), R? = 0.982 (20) 


Ga/Go = 5.292-4.226K;-0.321(S/Sp), R? = 0.991 (21) 


Eqs. (16)-(21) were used to calculate (Gq) and the obtained 
results were compared with the measured values of (Ga). The 
accuracy of estimating (Hg) was checked by calculating the MBE, 
RMSE, MPE, R°? and the t-test. 

A simple, physically based method proposed in [3,85,86] was 
used for estimating hourly diffuse and direct components from 
hourly total irradiance. For three different ranges of atmospheric 
transmissivity (K,=Gp/Gon), the resulting correlations are given by 
the following expression: 


0.995—0.081K; if Ke < 0.21. 
Dn/Gn = ¢ 0.724 + 2.738K,—8.32K? if 0.21 <K; <0.76. (22) 
0.180 if Kı > 0.76. 


Then, we can calculate the value of the hourly direct solar 
radiation as follows: 


Bn = Gh—Dh (23) 


4.3. Global solar energy models on inclined surface 


The knowledge of the solar total radiation incident on such a 
tilted surface is a prerequisite for the design of cost effective 
systems. The amount of elevation from the horizontal, the tilt 
angle, should be equal to the latitude angle of the location of the 
collector. This orientation is often selected for flat-plate collector 
installations, since it averages the installation peaks over the year. 
The surface axis tilted from the horizon by the latitude angle 
toward the south, is called an equatorial mounting. This is the 
ideal way of setting up the collectors. It is also the easiest and 
cheapest [3,10]. 

The total solar irradiation on a horizontal surface has been 
measured in many meteorological stations around the world, but 
there are only a few stations that measure the solar component on 
inclined surfaces. There are a number of models available to 
estimate solar irradiation on an inclined surface from the radiation 
on a horizontal surface, but these models require knowledge of the 
total irradiation, direct or diffused irradiation or reflected irradia- 
tion on a horizontal surface. An inclined surface solar irradiation 
model developed by Olmo et al. [3,10,44] requires only the 
horizontal surface solar irradiation, with incidence and solar 
zenith angles as input parameters. 

The Olmo et al. [3,44] model was developed to estimate the 
solar radiation on inclined surfaces using the data collected on 
horizontal surfaces. This model depends on the clearness index 
and avoids the direct and diffused solar radiation components. In 
the case of no ground reflections, the Olmo et al. model estimates 
the global irradiance (Gs) on an inclined surface from the 


corresponding solar radiation (G) on a horizontal surface by the 
following equation: 


Gp = Gy, (24) 


where £ is the surface inclination angle and ¥ is a function that 
converts the horizontal solar radiation to that incident on a tilted 
surface and is given as 


Wo = exp[-K(0"—-02)] (25) 


where @ and 6z (in radians) are the incidence and solar zenith 
angles, respectively, and K, is the hourly clearness index. 

Further, Olmo et al. [44] proposed a multiplying factor (F,) to 
take into account anisotropic reflections and it is given as 


Fe=1+psin7(0/2) (26) 


where p is the albedo of the underlying surface; this is the most 
commonly used expression for the radiation reflected from the 
ground. In this work a constant value for the albedo is used equal 
to 0.2. 

The Olmo et al. model for determining the global solar radia- 
tion on an inclined surface from that on a horizontal surface is 
then 


Gp = Gy,F. (27) 


The hourly total solar irradiance incident on a tilted surface 
(Grn) can be divided into three components: the beam component 
from direct irradiation of the tilted surface (Bra) and the ground 
reflected (Rr) and sky-diffuse (Drn) components: 


Grn = Bry + Dry + Rtn (28) 


The amount of direct radiation on a surface tilted S degrees 
from the horizontal and rotated a; degrees from the north-south 
axis can be calculated by multiplying the direct horizontal irradia- 
tion by the ratio of cos 6/cos Z, where @ is the solar incidence angle 
on a tilted plane and Z is the solar zenith angle. Also, the 
measuring station was located on a roof-top with very low 
reflectance, and the reflected component was very much lower 
than the direct and the diffuse components so an isotropic model 
can be used to compute the reflected component on the tilted 
surface. So, Eq. (28) can be written again as follows: 


Grn = By cos (0)/ cos (Z) + RaDn + Gpp[1- cos (S)]/2 (29) 


where By, Dn and Gn are the hourly direct, diffuse and total solar 
radiation on a horizontal surface, either measured directly or 
estimated from each other, Rg is the ratio of the hourly diffuse 
irradiation incident on a tilted surface to that on a horizontal 
surface, p is the ground reflectivity equal to 0.2 as assumed here 
because of the lack of a specific measurement, Z is defined from 
Eq. (6) and @ is calculated by the formula [86]: 


cos (Z) = sin (S) sin (Z) cos (aS—aT) + cos (S) cos (Z) (30) 


Several published meteorological data give the total solar 
radiation on horizontal surfaces; correlation procedures are 
required to obtain insolation values on tilted surfaces from 
horizontal radiation. Monthly average daily total radiation on a 
tilted surface (Gr) is normally estimated by individually consider- 
ing the direct beam (Gg), and diffuse (Gp) and reflected compo- 
nents (Gr) of the radiation on a tilted surface. Thus for a surface 
tilted at a slope angle from the horizontal, the incident total 
radiation is given by the relation: 


Gr = Gg + Gp + Gr (31) 


Several models have been proposed by various authors [3,6,30- 
35,87-89] to calculate solar radiation on tilted surfaces from the 
available data on a horizontal surface. The only difference among 
the models appears in the assessment of sky-diffuse component. 
Based on the assumptions made, the estimation models can be 
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classified into isotropic [90] and anisotropic [87,88] ones. The daily 
beam radiation received on an inclined surface can be expressed as 


Gg = (Gg—Ga)Rp (32) 


where Gg and Gg are the monthly mean daily global and diffuse 
radiation on a horizontal surface, and R, is the ratio of the average 
daily beam radiation on a tilted surface to that on a horizontal 
surface. The daily ground reflected radiation can be written as 


(1—cos p) 
2 


where £ is the tilt angle of the solar panel. Liu and Jordan [6,89] 
have suggested that R» can be estimated by assuming that it has 
the value which would be obtained if there was no atmosphere. 
For surfaces in the northern hemisphere, sloped towards the 
equator, the equation for R», is given as [89] 


Gr = Ggp (33) 


R cos (p—p) COS 5SiN øss + wss SiN (~—f) SiN ô (34) 
y= COS pcos SSiINwss + wss SiN osin ô 


where ass is the sunset hour angle for tilted surface for the mean 
day of the month, given by Eq. (3). For surfaces in the southern 
hemisphere, sloped towards the equator, the equation for Rp is 
given as follows: 


cos (p + P) cos ŝsin wss + wss Sin (p + p) sin 6 
COS PCOS SİN @ss + Mss SİN sin ô 


Rp = (35) 
and then, the total solar radiation on a tilted surface can be 
expressed as follows: 


(1—cos p) 
2 


As no information is available on ground albedo, p values are 
assumed to be 0.2; according to Eq. (36), we need the direct and 
diffuse components of global radiation for the estimation of global 
solar radiation on tilted surfaces [6]. 


Gr = (Gg—Ga)Rp + Ggp + GaRa (36) 


4.4. Diffuse solar energy models on inclined surface 


The models used to estimate the ratio of diffuse solar radiation 
on a tilted surface to that on a horizontal surface are classified as 
isotropic and anisotropic models. The isotropic models assume 
that the intensity of diffuse sky radiation is uniform over the sky 
dome. Hence, the diffuse radiation incident on a tilted surface 
depends on the fraction of the sky dome seen by it. The anisotropic 
models assume the anisotropy of the diffuse sky radiation in the 
circumsolar region (sky near the solar disk) plus an isotropically 
distributed diffuse component from the rest of the sky dome. The 
sky-diffuse solar radiation can be expressed as 


Ga = RaGa (37) 


where Rg is the ratio of the average daily diffuse solar radiation on 
a tilted surface, to that on a horizontal surface; the diffuse solar 
radiation models chosen in the present study are as follows 
[3,6,30]. 

The isotropic sky model [4,91-93] is the simplest model that 
assumes all diffuse solar radiation is uniformly distributed over 
the sky dome, i.e., it is independent of the azimuth and zenith 
angles. It approximates the completely overcast sky condition. The 
formula for the hourly sky diffuse solar radiation incident on an 
inclined plane is given by the product of the hourly diffuse solar 
radiation incident on a horizontal surface and the configuration 
factor from the surface to the sky, (1+cos #)/2. For the surface 
tilted by an angle p from the horizontal plane, the total solar 
irradiance can be written as follows: 


(1— cos £p) 


7 (38) 


Gar = Gay 


The isotropic model is given by the Badescu model (Ba) [94] as 
follows: 


_ 3+ cos (3p) 


Ra A (39) 
The Tian et al. model (Ti) [95] is given by the relation: 

Ra = 1—£/180 (40) 
The Koronakis model (Kr) [96] is given as follows: 

Ra =1/3[2 + cos(f)] (41) 
The Liu and Jordan (LJ) model [90] is given as follows: 

me —— (42) 


In addition to isotropic diffuse and circumsolar radiation, the 
Reindl model also accounts for horizon brightening [97,98] and 
employs the same definition of the anisotropic model. 

The Reindl et al. model (Re) [98] is given by the relation: 


Ra = (Gp/Go)Rp 


+|1-Go/Go)| SFI [1 + Vosne] (43) 


The Skartveit and Olseth model (SO) [99] is given as follows: 


Ra = (Gb /Go)Rp + 2 cos B + [1-(Gy/Go)—2] £ + 2 4 (44) 
where 
Q = {Max[0(0.3-2G; /Go)] (45) 


The Steven and Unswoth model (SU) [100] is given as follows: 
Ra = 0.51Rp 


+14 S050 _(4.74/1.26n)| sin p—p(x/180) cos p—x sin ?(/2)] 
(46) 
The Hay model (Ha) [91] is given by the relation: 
1 
Ra = (Gp/Go)Ry + [1—(Go/Ge)] [=F (47) 


The Klucher model, 1979 (K1) [101], is based on a study of clear 
sky conditions by Temps and Coulson [102]; their model was 
modified by Klucher, who incorporated conditions of cloudy skies. 
Klucher's formulation of the hourly sky diffuse solar radiation 
incident on an inclined surface is 


Gar = Gae{ [U T =R 


| [1 + Fy sin*(B/2)| [1 +F,cos76, sino,] } 
(48) 


where F; is the modulating function given by F, =(GaglGg)’, when 
the skies are completely overcast, F=0, Klucher's model reverts to 
the isotropic model. 

The Perez model (Perez et al. [103,104]; P8 and P9 respectively) 
[4] is more computationally intensive and represents a more 
detailed analysis of the isotropic diffuse, circumsolar and horizon 
brightening radiation by using empirically derived coefficients. 
The total irradiance on the tilted surface is given by the following 
equation: 


1+ cos f a 1—cos 
Gr = Gn bR + Gha |0 Fy)! P HF 2) 


5) ip Fasin 4 t Gio 5 
(49) 


here, F; and F, are the circumsolar and horizon brightness 
coefficients, respectively, and (a) and (b) are the terms that take 
the incidence angle of the sun on the considered slope into 
account. The terms (F), (F2), (a) and (b) are computed using the 
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following equations: 


nOZ noZ 
F= max 0, (fn +fi2d + Teo) & Fy =fo1 + food + g0 2 


(50) 


a= max(0°, cos @ and b = max(cos85, cos 02) (51) 


The coefficients f11, fiz, fis, for, f22, and fo3 were derived based on 
a Statistical analysis of empirical data for specific locations. Two 
different sets of coefficients were derived for this model [102,103]. 


5. Comparison techniques of modeling 


The relative ability of the different models to predict the solar 
radiation on horizontal and tilted surfaces was tested. The perfor- 
mance of the individual models was determined by utilizing 
statistical methods. There are numerous works in literature which 
deal with the assessment and comparison of daily solar radiation 
estimation models. The most popular statistical methods are the 
MBE (mean bias error) and the RMSE (root mean square error). In 
this study, to evaluate the accuracy of the estimated data, from the 
models described above, the following statistical estimators were 
used: MBE, RMSE, MPE (mean percentage error) and the correla- 
tion coefficient (R°), to test the linear relationship between 
predicted and measured values. For higher modeling accuracy, 
these estimators should be closer to zero, and the correlation 
coefficient, (R?), should approach to 1. The NSE (Nash-Sutcliffe 
equation) was also selected as an evaluation criterion. A model is 
more efficient when NSE is closer to 1. However, these estimated 
errors provide reasonable criteria to compare models but do not 
objectively indicate whether the estimates from a model are 
statistically significant. The t-test statistic allows models to be 
compared and at the same time it indicates whether or not a 
model's estimate is statistically significant at a particular con- 
fidence level. So, the t-test was carried out on the models to 
determine the statistical significance of the predicted values [3,10]. 


5.1. Mean Bias Error (MBE) 


To evaluate the accuracy of the prediction data from the models 
described above, this test provides information on the long-term 
performance of a model. A low MBE value is desired. A negative 
value gives the average amount of underestimation in the calcu- 
lated value. So, one drawback of MBE is that overestimation of an 
individual observation may cancel underestimation in a separate 
observation. We can obtain the values of MBE as follows: 


1 
MBE = n Di (Gi-catc.—Gimeas.) (52) 
and the equation of mean percentage error MPE% is expressed by 
1 
MPE% = n Xi [(Gicaic — Gimeas.)/Gimeas. 100 (53) 


The subscript i refers to the ith value of the daily solar 
irradiation and n is the number of the daily solar irradiation data. 
The subscripts “calc.” and “meas.” refer to the calculated and 
measured daily solar irradiation values, respectively. A percentage 
error between — 10% and +10% is considered acceptable [105]. 


5.2. Root Mean Square Error (RMSE) 


The value of RMSE is always positive, representing zero in the 
ideal case. The normalized RMSE gives information on the short- 
term performance of the correlations by allowing a term-by-term 
comparison of the actual deviation between the predicted and 
measured values. The smaller the value is, the better the model's 


performance is, and the equation of RMSE is as follows [106]: 


1 


1/2 
RMSE = i E? Gica —Gimeas)]| (54) 


5.3. The t-test statistic (t) 


The tests for mean values, the random variable t with n—1 
degrees of freedom may be written here as follows [107]: 


t = [(n—1)(MBE)? /(RMSE)* —(MBE)*]'/2 (55) 


The smaller the values of t-statistic the better the performance 
of modeling. 


5.4. The correlation coefficient (R°) 


In statistics literature, it is the proportion of variability in a data 
set that is accounted for by a statistical model, where the 
variability is measured quantitatively as the sum of square devia- 
tions. Most often it is defined notationally as [108] 


R = YP (Xi-YH21/LP_ [(Xi-Yi)2] (56) 
This can also be expressed as 
R = 1—(YP_ [Xi-YH2]/LP_ {[Xi-Yi)2)} (0 <R° <1) (57) 


herein, Xi and Yi are the measurements and model estimates, 
respectively. A high value of R? is desirable as this shows a lower 
unexplained variation. R? is a statistic that gives some information 
about the goodness-of-fit of a model. In regression, the R? 
coefficient of determination is a statistical measure of how well 
the regression line approximates the real data points. An R? of 
1.0 indicates that the regression line perfectly fits the data, which 
is never valid in any solar radiation estimation model. 


6. Results and discussion 


The regression equations between (G/G,) and meteorological 
variables along with the values of MBE, RMSE, MPE, R? and the t- 
test statistics are summarized in Table (1). The empirical correla- 
tions for estimation of total solar irradiation on horizontal surface 
in the form of Eqs. (7)-(15) are proposed for Cairo using the 
meteorological data during the time period 1990-2010. From the 
analysis of the measured and calculated values of the total solar 
irradiation (G), it is clear that, the values of correlation coefficients 
(R?) are higher than 0.95 and the values of the RMSE are found in 
the range 3.13-6.34, thus indicating a good agreement between 
measured and calculated values of the total solar radiation (G). The 
negative values of the MPE show that, Eqs. (8), (9), (12), and (13) 
slightly overestimate the values of the total solar irradiation (G), 


Table 1 
The relation between (G/G,), (W/m?) and metrological variable at Cairo, Egypt, 
during the time period (1990-2010) by using the empirical correlations models. 


Model no. Regression coefficients MBE RMSE MPE% R? t-test 


a b c 
Eq. (7) 0.546 0.356 - —3.43 5.17 6.55 0.978 5.45 
Eq. (8) 0.578 0.389 0.613 —4.43 425 -5.34 0.986 4.67 
Eq. (9) 0.243 0.456 0.676 —2.34 498 -2.92 0.992 4.32 
Eq. (10) 0.278 0.675 0.587 —3.87 3.67 2.65 0.965 3.45 
Eq. (11) 0.353 —0.439 0.676 2.76 3.16 5.23 0.965 3.67 
Eq. (12) 0.316 0.543 0.434 6.56 554 —3.76 0.984 6.34 
Eq. (13) —0.487 0.436 —0.687 —4.89 5.89 -2.36 0.973 7.65 
Eq. (14) —0.346 0.489 0.422 —3.21 3.11 6.54 0.962 2.87 
Eq. (15) 0.523 0.354 0.632 —2.43 435 7.21 0.964 5.65 


Table 2 
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Differences between measured (Gam) and calculated (Ga,c) values (W/m?) with the 


metrological variable at Cairo, Egypt, during the time period (1990-2010). 


Model no. 


Table 3 


Gam 


1347 
1454 
1421 
1362 
1348 
1351 


Gaye 


1387 
1496 
1374 
1329 
1378 
1331 


MBE 


— 1.87 
3.65 
2.42 
1.89 

—2.75 

—2.47 


RMSE 


2.38 
2.91 
2.24 
2.47 
1:21 
2.68 


MPE% 


—1.89 
2.62 
1.91 
1:25 
1.39 

—2.74 


Statistical results analysis for Olmo et al., in the present study. 


Month 


January 
February 
March 
April 

May 

June 

July 
August 
September 
October 
November 
December 


Table 4 


MBE 


29.45 
24.82 
19.32 
18.71 
(7.34 
9.34 
8.84 
6.78 
8.12 
13.24 
19.76 
32.57 


RMSE 


42.24 
37.57 
31.34 
30.24 
28.78 
22.76 
19.28 
22.61 
20.87 
26.47 
33.91 
45.18 


MPE% 


9.25 
6.97 
6.34 
5.78 
4.87 
2.56 
3.89 
4.46 
5.78 
6.55 
7.46 
8.78 


R? 


0.934 
0.958 
0.991 
0.982 
0.974 
0.991 
0.984 
0.986 
0.969 
0.986 
0.942 
0.926 


t-test 


2.65 
3.24 
1.78 
1.16 
1.14 
2.16 


t-test 


4.89 
4.24 
3.45 
4.74 
4.25 
3.14 
2.86 
3.54 
4.37 
6.45 
5.24 
7.47 


but, Eqs. (7), (10), (11), (14) and (15) slightly underestimate the 
values of the global solar irradiation (G). These results indicate a 
good agreement between the average hourly daily data of the total 
solar radiation and the other meteorological parameters. Also from 
Table 1, it is seen that, the values of (t-test) change from one model 
to another model according to the models from Eqs. (7)-(15). Thus 
the model which gives the smallest values of the t-test is then 
considered as the best model for estimating the total solar 
irradiation at the selected site with an acceptable error. This 
means that the models of Eqs. (10), (11) and (14) are a good 
estimate for the total solar irradiation in the selected location 
during the time period in the present work. 

Table 2 shows the differences between the measured (Gam) and 
calculated (Gac) values of the diffuse solar radiation along with the 
values of mean base error (MBE), root mean square error (RMSE), 
mean percentage error (MPE), and t-test statistics. From this table, it 
is clear that, the low values of the RMSE for all models indicate a 
good agreement between measured and calculated values of diffuse 
solar radiation (Gg). The negative values of MPE indicate that the 
proposed correlations slightly overestimate (Gg). For all models, the 
absolute values of the MPE indicate very good agreement between 
measured and calculated values of the diffuse solar fraction (Gg/G) or 
the diffuse solar transmittance (Ga/Go) and clearness index K, 
relative number of sunshine hours (S/S,) and their combination. 
Also from Table 2, the t-test of the model in Eqs. (19) and (20) is 
given the smallest value, and so it is considered as the best model for 
estimating the diffuse solar radiation at the selected site with an 


Average daily values of the mean base error, MBE (%), for global solar irradiation received on an inclined surface at different slopes in the time period 1990-2010. 


Slope Jan. Feb. 
Isotropic model MBE 
0 — 0.184 — 0.356 
15 — 4.48 —3.87 
30° —7.68 —5.77 
45° -911 — 11.42 
60° — 10.35 — 12.54 
75 —9.43 — 11.67 
90° —8.12 —11.76 
Hay's model MBE 
0 — 0.176 — 0.223 
15° —3.59 —2.87 
30° — 4.96 —4.23 
45° — 4.34 —3.87 
60° —7.76 —6.95 
75° —7.11 —7.54 
90° —4.78 —5.84 
Klucher's model MBE 
0 — 0.187 —0.235 
15 — 2.65 —2.69 
30° —4.76 —3.87 
45° —3.89 —2.94 
60° —5.45 —4.87 
75° —4.78 —5.57 
90° —3.76 —4.65 
Perez model MBE 
0 — 0.154 —0.275 
15 —2.45 —2.65 
30° —2.24 — 2.89 
45° —3.87 —3.47 
60° —1.98 — 2.76 
75° —2.32 —2.57 
90° — 4.46 —3.34 
Temps and Coulson's models MBE 
0 2.43 2.68 
15 3.76 3.32 
30° —1.78 2.32 
45° —2.56 2.45 
60° —5.76 —5.16 
75° — 4.94 —3.89 
90° —5.73 —4.23 


Mar. 


—0.164 


-3.65 
—4.34 
—5.49 


—5.64 


—0.287 


—0.234 
—2.49 
—1.98 
—2.41 
—2.87 
—2.89 
—3.45 


3.73 
2.56 
2.87 
3.43 
—4.76 
— —2,98 
—3.54 


Apr. 


—0.198 
—1.56 
—3.32 
—5.41 
—5.89 
—6.56 
—5.76 


—0.186 
—1.76 
—2.87 
—3.45 
—5.93 
— 4.87 
—3.56 


—0.154 
— 1.0.23 
—2.65 
—1.76 
—2.92 
—3.76 
—2.76 


—0.147 
—1.86 
—2.37 
—2.74 
—1.87 
—3.15 
—2.67 


2.87 
3.32 
4.87 
3.76 
2.95 
1.87 
3:32 


May Jun. July Aug. Sept. Oct. Nov. Dec. 
0.256 0.065 0.184 0.148 0.224 0.176 0.132 0.189 
1.56 0.76 0.69 1.54 2.78 2.34 3.79 3.78 
2.43 1.21 1.54 2.61 4.76 4.43 8.55 7.56 
3.43 2.34 0.76 3.76 5.32 6.78 9.23 6.23 
5.34 4.23 2.78 3.55 6.67 7.92 7.65 8.89 

—6.1 3.15 1.97 4.23 7.84 8.43 6.88 9.56 

—3.69 1.87 1.54 2.54 3.35 5.54 7.43 8.19 
0.245 0.132 0.129 0.165 0.123 0.143 0.154 0.134 
0.69 0.97 0.76 1.43 1.47 2.67 2.87 2.34 
2.21 0.78 0.54 1.78 3.87 2.34 412 3.34 
2.78 1.32 1.87 2.54 3.23 3.76 5.54 4.21 
4.54 2.34 1.89 2.68 3.56 5.4 6.43 6.85 
3.32 1.87 1.48 2.87 3.23 4.23 5.67 5.93 
1.47 1.32 0.115 1,91 2.56 3.41 3.65 4.47 
0.286 0.113 0.149 0.167 0.146 0.178 0.176 0.134 
0.79 0.89 0.87 1.54 1.43 2.67 2.32 2.86 

—1.45 0.134 0.46 1.49 2.67 2.87 3.76 3.21 

— 1.78 1.56 1.78 2.76 2.32 3.79 2.95 3.85 

—1.59 2.76 1.33 2.32 3.54 4.76 5.76 5.21 

—2.76 1.86 1.74 2.76 2.76 3:91 4.86 4.56 

—1.84 1:32 1.65 2.89 2.45 3.56 4.43 3.86 
0.267 0.124 0.146 0.184 0.124 0.158 0.184 0.132 
0.98 0.86 1.27 1.46 1.34 2.79 2.94 2.24 
1.64 0.49 0.85 1.68 1.79 1.67 2.54 2.41 

— 1.87 —1.48 0.59 2.24 1.43 1.24 2.24 2.69 

—2.45 1.68 1.29 1.67 2.78 2.84 2.79 2.23 

— 2.87 2.14 1.84 1.74 2.14 2.39 3.45 1.98 

-1.23 1.85 1.47 2.39 3.29 3.47 4.21 5.65 
3.48 2.65 2.87 3.76 3.25 4.43 4.35 5.76 
2.34 1.89 3.54 3.97 2.56 3.76 4.87 4.34 
5.65 6.76 4.78 5.56 3.22 —3.87 —2.43 2.78 
4.87 5.76 5.23 —4.87 —4.89 2.32 1.95 3.34 
6.56 8.43 4.86 5.21 3.76 2.18 2.86 2.94 
2.44 4.32 3.54 3.34 1.87 5.76 3.76 1.65 
7.54 8.36 — 4.93 —4.87 4.32 4.98 1.24 — 2.78 
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acceptable error. These results are in good agreement with the 
previous work performed for other studies [3,57]. 

The results of the statistical analysis of the relative ability of the 
Olmo et al. model to determine the solar global irradiation on the 
inclined surface are presented in (Table 3). The Olmo et al. 
{3,10,44] model was applied to the database corresponding to 
the horizontal solar total irradiation to determine values for a 
south facing surface, tilted at latitude angle for all sky conditions. 
The hourly daily solar irradiation data have been used in this study 
and the statistical coefficients have been computed on the basis of 
the experimental data. From Table 3, we notice that, the mean 
percentage error (MPE) is in the range between 2.56% and 9.25%, 
but the values of root mean square error (RMSE) varies from 19.28 
to 45.18. Also from this table, we notice that all values of MBE and 
RMSE in all months are nearest to these values, with the exception 
of winter months. Normally, these months are subject to heavy 
rainfall and less solar radiation. The high values of MPE and RMSE 
and low value of correlation coefficient for these months can be 
justified. The model provides a good estimation tool for the other 
months. In general, considering the statistics as a whole, the total 
solar irradiation data estimated by the Olmo et al. model are in 
good agreement with the measured values. Therefore, the Olmo 
et al. model is recommended to estimate the total solar radiation 
on an inclined surface in this study, due to its accuracy, input 
requirements and simplicity. 


Table 5 
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The average hourly daily values of the mean base error, MBE 
(%), for global solar irradiation received on the inclined surface at 
different slopes during the time period 1990-2010 are summar- 
ized in Table 4. From this table, we notice that, five solar irradiance 
models are considered in this work: the isotropic model, and the 
Hay, Klucher's, the Perez, and Temps and Coulson's which are the 
anisotropic models. For each model, the measured values of diffuse 
solar radiation and horizontal values of total solar radiation were 
used to calculate the solar radiation on the surface tilted at 
different slopes which varies from 0° to 90° above the horizon. 
The results were compared with the solar irradiances monitored 
and presented in terms of usual statistics: the mean base error 
(MBE) and the root mean square error (RMSE). Also from Table 4, it 
is seen that, the values of MBE vary from —0.065 to —12.54, 
— 0.123 to —7.76, — 0.113 to —5.76, — 0.124 to —4.46 and 1.24 to 
8.43 for the isotropic, Hay, Klucher, Perez and Temps and Coulson's 
models respectively. The values of MBE results show that, the 
isotropic, Perez's, Hay's and Klucher's models substantially under- 
predict the irradiance incident on an inclined surface, and the 
Temps and Coulson model considerably overpredicts irradiance 
incident on an inclined surface on an overall basis. 

Table 5 shows the results for south facing surfaces of the root 
mean square error (RMSE) for different slopes in different models 
in the present work. From this table, we observed that, the values 
of root mean square error vary from 0.31 to 31, 0.29 to 21, 0.25 to 


Average daily values of the root mean square error, RMSE (%), for global solar irradiation received on an inclined surface at different slopes in the time period 1990-2010. 


Slope Jan. Feb. Mar. Apr. May 
Isotropic model RMSE 
0 0.79 0.49 0.67 0.84 0.76 
15 14 15 10 9 7 
30° 18 19 13.7 11 8 
45° 22 27 25 16 12 
60° 17 21 27 19 19 
75 19 25 22 26 26 
90 23 31 24 21 17 
Hay's model RMSE 
0 0.76 0.49 0.67 0.34 0.72 
15 9 12 7 6 4 
30° 11 15 9 8 65 
45° 9 8 12 5 6 
60° 8 14 8.5 9 9 
75 13 16 11 12 7 
90° 16 19 13 9 11.5 
Klucher's model RMSE 
0 0.49 0.56 0.76 0.37 0.54 
15 6 71 3.9 4.6 3.6 
30° 9.3 10.5 6.8 5.4 7 
45 73 7.6 9 4 5.2 
60° 4.6 9.6 72 7 7 
75° 8 12 11 10.5 9.3 
90° 13 15.6 8.9 9 11.7 
Perez model RMSE 
0 0.44 0.37 0.25 0.19 0.13 
15 1.65 2.45 2.87 3.65 2.54 
30 1.89 2.46 2.47 2.86 3.4 
45 3 3.76 2.54 2.43 1.78 
60° 5 7.89 11 6 5.76 
75 8.56 10.4 9.4 9.2 7.2 
90° 11 9.3 8 7.2 5.8 
Temps and Coulson's models RMSE 
0 3.47 5 4.5 6 7 
15 5.4 6.3 8.3 9 10.5 
30° 7 75 10.4 11 12 
45 8.3 9 9.2 12 14 
60° 10.5 12 15 14 16 
75 11.7 13 16 17 19 
90° 14 15.7 18 20 22 


Jun. 


July Aug. Sept. Oct. Nov. Dec. 
0.31 0.54 0.87 0.67 0.47 0.54 
5 6 6.7 5 7 11.3 
6 7.8 9.2 11 13 19 
12 11 5.6 17 24 31 
17 14 10 9 17 26 
16 18 14 16 25 27 
19 14 24 22 30 24 
0.29 0.39 0.56 0.64 0.32 0.48 
3.4 42 48 5.3 5.2 6.8 
2.3 3.3 6 9 8.4 11.8 
5 6.2 43 75 10 17.2 
3.7 7.8 6.4 6.7 7.8 14 

7 10.6 10 12.6 10.4 13.5 
8.7 12 13 15.2 12 12.2 
0.25 0.32 0.52 0.42 0.37 0.63 
2.3 3.8 4 4.7 5.8 6.9 
17 43 6 7A 7A 9 
49 5.6 5.7 8.6 10 12 

6 8 8.2 6 75 15 
8.5 10.7 12.5 10.3 9.6 13.2 
9.7 12.4 15.4 12 11.5 11.4 
0.08 0.55 0.37 0.34 0.24 0.32 
1.57 2.65 2.11 3.76 223 2.87 
1.34 3.49 2.78 2.23 3.65 4.11 
1.54 3.23 3.23 3.76 3.45 5.23 
46 6.4 7.45 8.45 732 8.56 
6 5.8 8.46 7.34 9.43 9.23 
8.4 7.45 8.56 10.4 11.6 11.4 
10 7 55 46 5 75 
12 10 7 7 8.6 9 

13 12 9 9 11 11 
15 14 11 10 13 14 
18 21 18 15 17 17 
26 26 25 23 19 22 
30 31 28 30 28 25 


S.A. Khalil, A.M. Shaffie / Renewable and Sustainable Energy Reviews 27 (2013) 853-863 861 


Table 6 


The statistical results of the models for the south-facing surface at Cairo, Egypt, during the time period from 1990 to 2010. 


Model Abbreviation Year a 

Liu and Jordan LJ (ISO) 1962 0.154 
Temps and Coulson TC (ANI) 1977 0.294 
Hay Ha (ANI) 1979 0.345 
Stevenand Unsworth SU (ANI) 1980 0.534 
Skartveit and Olseth SO (ANI) 1986 0.674 
Koronakis Kr (ISO) 1986 0.165 
Perez et al. P8 (ANI) 1986 0.281 
Reindl et al. Re (ANI) 1990 0.347 
Perez et al. P9 (ANI) 1990 0.449 
Tian et al. Ti (ISO) 2001 0.487 
Badescu Ba (ISO) 2002 0.565 


ISO means isotropic and ANI means anisotropic. 


Table 7 


The statistical results of the models for the west-facing surface at Cairo, Egypt, during the time period from 1990 to 2010. 


Model Abbreviation Year a 

Liu and Jordan LJ (ISO) 1962 0.423 
Temps and Coulson TC (ANI) 1977 0.312 
Hay Ha (ANI) 1979 0.396 
Steven and Unsworth SU (ANI) 1980 0.546 
Skartveit and Olseth SO (ANI) 1986 0.497 
Koronakis Kr (ISO) 1986 0.411 
Perez et al. P8 (ANI) 1986 0.511 
Reindl et al. Re (ANI) 1990 0.478 
Perez et al. P9 (ANI) 1990 0.396 
Tian et al. Ti (ISO) 2001 0.432 
Badescu Ba (ISO) 2002 0.532 


ISO means isotropic and ANI means anisotropic. 


15.6, 0.06 to 11 and 3.47 to 30 for the isotropic, Hay, Klucher, Perez 
and Temps and Coulson's models respectively. From these values 
we conclude that, the RMSE values for the fifth model increase as 
the slope of the collector increases, but remain in a domain of 
error for which these relations can be applied with good accuracy. 
Inspecting the results, it is apparent that the models agree quite 
well with each other during the summer months. They deviate 
from each other in the winter months, when the effect of the 
difference in the diffuse solar radiation parameterization is at its 
maximum. The RMSE results indicate that the anisotropic models 
(Hay, Klucher and Perez) show similar performance on an overall 
basis, but isotropic model and Temps and Coulson model exhibit 
much larger error. In general we confirm that, the observation of 
the Perez and Klucher models describes the irradiance on inclined 
plane more accurately than other models. These results in the 
present work are in good agreement with the other works which 
are found in [3,108]. 

Eleven models have been used to evaluate the statistical 
parameters: mean base error (MBE), root mean square error 
(RMSE), mean percentage error (MPE), correlation coefficient (R°) 
and t-test statistics for south-facing and west-facing surfaces in 
this paper which can be referred to in Tables 6 and 7. The models 
used were the isotropic models of Badescu (Ba) [94], Tian et al. (Ti) 
[95], Koronakis (Kr) [96] and Liu and Jordan (LJ) [90], and the 
anisotropic models of Perez et al. (P9) [104], Reindl et al. (Re) [98], 
Perez et al. (P8) [103], Skartveit and Olseth (SO) [99], Steven and 
Unsworth (SU) [100], Hay (Ha) [87], Klucher (K1) [101], and Temps 
and Coulson (TC) [102]. 

The evaluation was carried out on a semi-hourly basis. The total 
solar irradiation component on the tilted surface was determined 
from measured horizontal data using different models and com- 
pared with the measured tilted data of the same period in the 
present work. Tables 6 and 7 report a summary of the statistical 
results of the models for south facing and west facing surfaces in 


b MBE RMSE MPE R? t-test 
0.425 12.6 28.7 —8.4 0.957 7.57 
0.656 —16.4 33.3 11.8 0.974 8.41 

0.511 —30.4 9.3 14.6 0.972 6.11 

0.632 25.8 42.5 16.7 0.981 11.74 
0.341 —13.6 19.4 -11.7 0.974 6.47 
0.753 -27.9 25.9 13.7 0.962 8.11 

0.524 18.4 17.9 16.9 0.958 8.87 
0.524 19.6 8.7 12.4 0.991 4.64 
0.319 —15.8 19.5 -9.3 0.974 5,55 
0.384 18.7 24.2 12.8 0.985 7.89 
0.376 —22.4 29.5 11.6 0.942 8.47 
b MBE RMSE MPE RP t-test 
0.342 —12.4 34.2 -13.5 0.974 8.45 
0.487 —15.2 46.3 16.8 0.949 11.54 
0.498 —35.6 19.9 19.2 0.987 7.92 
0.343 11.3 47.3 10.8 0,953 13.76 
0.367 14.7 22.5 —14.9 0.965 6.34 
0.398 —25.7 31.2 17.6 0.971 8.94 
0.332 16.9 25.8 19.4 0.978 8.23 
0.322 —19.3 19.7 13.8 0.983 6.64 
0.456 — 12.6 12.4 12.4 0.992 4.76 
0.325 11.2 30.8 15.8 0.985 6.54 
0.387 —22.4 33.6 19.2 0.974 10.76 


the present study respectively. It is seen from Table 6 that the 
absolute relative values of the root mean square error (RMSE), for 
the south facing surface, range from 8.7 to 42.5 for the Hay (Ha) 
model and the Stevenand Unsworth (SU) model respectively. 
These results are in good agreement with the values of t-test 
statistics for the same models and the correlation coefficient is 
clear with the higher value too for self-models. For west-facing 
surface shown in Table 7, the values of root mean square error 
range from 12.4 to 47.3 for the Perez et al. (P9) and Temps and 
Coulson (TC) models respectively. From Tables 6 and 7, we can 
conclude that the Hay (Ha), Skartveit and Olseth (SO) and Perez 
et al. (P9) models give the most accurate predictions for the south- 
facing surface, and Hay (Ha) and Perez et al. (P9) models are 
perform better estimation for the west-facing surface. 


7. Conclusion 


The regression equations between (G/G,) and meteorological 
variables along with the values of MBE, RMSE, MPE, R? and the t- 
test statics are summarized in the present work. Empirical 
correlations for the estimation of total solar radiation on horizon- 
tal surface in the form of Eqs. (7)-(15) are proposed for Cairo using 
the meteorological data during the time period 1990-2010. From 
the analysis of the measured and calculated values of the total 
solar radiation (G), it is clear that, the values of correlation 
coefficients (R?) are higher than 0.95 and the values of the RMSE 
are found in the range 3.13-6.34, thus indicating a good agreement 
between measured and calculated values of the total solar radia- 
tion (G). The models of Eqs. (10), (11) and (14) give a good estimate 
of the total solar radiation in the selected location during the time 
period in the present work. 

The values of the RMSE for all models in Table 2 indicate a good 
agreement between measured and calculated values of diffuse 
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solar radiation (Ga). The negative values of MPE indicate that the 
proposed correlations slightly overestimate (Gg). For all models, 
the absolute values of the MPE indicate very good agreement 
between measured and calculated values of the diffuse solar 
fraction (Gy/G) or the diffuse solar transmittance (Ga/Go) and 
clearness index K, relative number of sunshine hours (S/S,) and 
their combination. 

The values of MBE results show that the isotropic, Perez's, Hay's 
and Klucher's models substantially underpredict the irradiance 
incident on an inclined surface, and the Temps and Coulson model 
considerably overpredicts irradiance incident on an inclined sur- 
face on an overall basis, and the RMSE results indicate that the 
anisotropic models (Hay, Klucher and Perez) show similar perfor- 
mance on an overall basis, but isotropic model and the Temps and 
Coulson model exhibit much larger error. In general we confirm 
that, the observation of Perez's and the Klucher models describes 
the irradiance on an inclined plane more accurately than other 
models. These results in the present work are in good agreement 
with the other works in [3,108]. 

The Hay (Ha), Skartveit and Olseth (SO) and Perez et al. (P9) 
models give the most accurate predictions for the south-facing 
surface, and Hay (Ha) and Perez et al. (P9) models perform better 
estimation for the west-facing surface. 
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